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Abstract 


The Multi-channel Photometric Survey Telescope (Mephisto) is a real-time, three-color photometric system 
designed to capture the color evolution of stars and transients accurately. This telescope system can be crucial in 
cosmological distance measurements of low-redshift (low-z, z <0.1) Type Ia supernovae (SNe Ia). To optimize the 
capabilities of this instrument, we perform a comprehensive simulation study before its official operation is 
scheduled to start. By considering the impact of atmospheric extinction, weather conditions, and the lunar phase at 
the observing site involving the instrumental features, we simulate light curves of SNe Ia obtained by Mephisto. 
The best strategy in the case of SN Ia cosmology is to take the image at an exposure time of 130 s with a cadence of 
3 days. In this condition, Mephisto can obtain hundreds of high-quality SNe Ia to achieve a distance measurement 
better than 4.5%. Given the on-time spectral classification and monitoring of the Lijiang 2.4m Telescope at the 
same observatory, Mephisto, in the whole operation, can significantly enrich the well-calibrated sample of 
supernovae at low-z and improve the calibration accuracy of high-z SNe Ia. 
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1. Introduction 


As the most reliable indicators of cosmological distance so 
far, Type Ia supernovae (SNe Ia) play a crucial role in 
unveiling the accelerated expansion of the universe (Riess et al. 
1998; Perlmutter et al. 1999). Additionally, SNe Ia have greatly 
contributed to the discovery of the Hubble tension (Dainotti 
et al. 2021). Despite the proposal of dark energy as an 
explanation for the accelerated expansion of the universe, our 
understanding of its fundamental properties remains limited. As 
the most straightforward observational tool for constraining the 
parameters associated with dark energy, SNe Ia achieve a high 
level of accuracy when employed for distance measurements 
(Weinberg et al. 2013). However, the recent accuracy of SNe Ia 
(e.g., 7%, Abbott et al. 2019) falls short in constraining 
cosmological parameters and determining the nature of dark 
energy. 

The crucial factor in constraining the range of cosmological 
parameter values lies in accurately ascertaining the distances of 
SNe Ia at high-z (z>0.5). The application of SNe Ia in 
cosmological distance measurement depends on some empiri- 
cal relations between the absolute peak brightness and shape of 
light curves, for example, the width-luminosity relation (WLR, 
usually known as the Phillips relation; Phillips 1993). The 


original form of the Phillips relation is 
M (A) max =at bAm (A), (1) 


where M (\)max represents the maximum absolute magnitude at 
a particular band (usually using the B-band), and Amy; is the 
decline rate from the peak luminosity to 15 days later. The 
coefficients a and b are measured by dozens of well-observed 
low-z SNe Ia. Comparing with the low-z sample, one can 
derive the intrinsic brightness of SNe Ia at a corresponding 
band for a given redshift and thus measure its luminosity 
distance. For example, at z = 0.01, the rest-frame 4310 A flux is 
observed in the B-band, while at z= 1.01 this same flux is 
roughly observed in the /-band. To accurately determine the 
distances of SNe Ia at various redshifts, it is imperative to 
possess knowledge regarding the wavelength-dependent evol- 
ution of WLR. This calibration predicament in SN Ia 
cosmology is the primary source of systematic error in 
constructing the SNe Ia Hubble diagram. Consequently, 
enhancing calibration precision at low-z becomes indispensable 
for refining measurements of cosmological parameters using 
SNe Ia. Moreover, low-z supernova samples not only help to 
reduce the uncertainties associated with high-z SN Ia measure- 
ments but also provide constraints to the Hubble constant 
(Riess et al. 2009, 2011, 2016, 2022), which can alleviate the 
Hubble tension. 
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For low-z samples, the current largest source of systematic 
error is the calibration. The existing low-z SNe Ia come from 
different telescopes, with some using Bessel-like UBVRI filters 
and some utilizing SDSS-like ugriz filters. Differing filter 
systems present a significant challenge for calibration at the 
level of a few percent. Even with the same filter system, subtle 
differences can exist in the transmissions. The most direct 
approach to reducing this error is to obtain additional samples 
using a consistent telescope system. Extinction is another factor 
of systematic uncertainty. For the purpose of calculating 
extinction, accurate color measurement is essential. Addition- 
ally, variations in ignition positions (central versus off-center) 
and explosion mechanisms (deflagration versus detonation) 
contribute to the diversity of SN Ia observations. 

Furthermore, WLR is not only a function of wavelength but 
also depends on the spectral parameters of SNe Ia. Besides the 
outliers that cannot be standardized by any kind of WLR, 
normal SNe Ia can be further divided into some subclasses 
depending on the spectral features around the peak brightness, 
e.g., the equivalent width (EW, Branch et al. 2006), velocity 
(Wang et al. 2009), and temporal velocity gradient (Benetti 
et al. 2005) of Sill 46355. The fine classification is an efficient 
way to improve the accuracy of SNe Ia in distance 
measurement. For example, Wang et al. (2009) divided normal 
SNe Ia into two subclasses depending on the velocity of Sill 
6355 around the maximum brightness and improved the 
accuracy of the Hubble diagram from 0.178 to 0.125 mag. 
Blondin et al. (2012) use different coefficients to fit each 
subclasses of various SNe Ia depending on the EW of Sill 
lines. Thus, to create a consistent sample, it is crucial to have 
three spectra for every well-observed SN Ia that covers the 
spectral evolution from a few days before to a few days after 
the B-band maximum. 

In this paper, we investigate the possibility of the Multi- 
channel Photometric Survey Telescope (Mephisto, Yuan et al. 
2020) for low-z calibration of SNe Ia. Mephisto is a 1.6 m 
wide-field high-precision telescope with a 2 deg? field of view. 
It is led by Yunnan University and co-built by Yunnan 
Observatories and Nanjing Institute of Astronomical Optics & 
Technology. This telescope can simultaneously image the sky 
in three bands (ugi or vrz) and obtain real-time color 
information. It boasts high color calibration accuracy, which 
can effectively reduce errors introduced by extinction and 
reddening. Moreover, the available samples around z= 0.1 so 
far are low (Abbott et al. 2019). As a 1.6 m-aperture wide-field 
survey telescope, it has the ability to obtain a great number of 
high-quality samples. Based on multi-color photometry, it can 
trigger the Lijiang 2.4m Telescope(LJT, Fan et al. 2015) to 
perform follow-up spectral observation. We can get the spectral 
evolution to help us identify the subclass. Hence, we can 
remove SNe Ia with larger deviations in the analysis, and thus 
reduce the uncertainty in the Hubble diagram. 
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To improve the observation efficiency of Mephisto, we 
conduct observation simulation and cosmological fitting to 
produce an appropriate observation strategy. We use SNCosmo 
(Barbary et al. 2022) in this study to simulate the light curves of 
SNe Ia for Mephisto. Section 2 provides a concise introduction 
to the fundamental information of Mephisto and the Lijiang 
Observatory. Section 3 delves into the simulation of SN Ia light 
curves and the photometric accuracy under different exposure 
durations. The cosmological parameter estimation method is 
detailed in Section 4. Section 5 is the presentation of the 
results, and there we have provided appropriate observation 
strategies. Section 6 presents the results of the estimation and 
provides suggestions for future observations. 


2. Mephisto and Lijiang Observatory 


Mephisto innovatively uses a Ritchey-Chrétien (RC) system 
with correctors and film-coated cubic prisms, so that the focal 
plane of the telescope is under the main mirror and there is a 
large distance between the focal plane and the main mirror. 
Therefore, the color separation system and multiple cameras 
can be placed in the rear optical path while ensuring a high- 
quality image in the entire field of view. The telescope is 
equipped with three high-quality CCD cameras, which can 
capture real-time color images of the same patch of sky in three 
different bands (ugi or vrz), providing valuable color 
information. 

The current total efficiency of Mephisto is depicted in 
Figure |. The filters used by Mephisto are not the most popular 
system among observers, e.g., UBVRI of Johnson—Cousins or 
ugriz of SDSS. The specific choice of these filters depends 
largely on their scientific objectives. Mephisto’s survey mode 
includes both W-survey and D/H/M-survey modes. The W 
mode focuses on 2*20 second exposures of observable areas in 
the northern hemisphere (6 > 15°, 26,000 square degrees). The 
D/H/M modes, on the other hand, sample at daily, hourly, and 
minute cadence, respectively, for time-domain astronomy 
observations. The observation strategy adopted in this paper 
is based on D mode. 

Mephisto is situated at the Lijiang Observatory administered 
by Yunnan Observatories, located in Gaomeigu, Lijiang City, 
in Yunnan Province, China, at an altitude of 3200 m. The site 
benefits from an effective observation time of more than 2000 
hr per year. The average seeing is about 1/17, with 25% of the 
time experiencing seeing below 0”91 and 75% of the time 
experiencing seeing below 1/33 (Xin et al. 2020). Overall, the 
Lijiang Observatory is an outstanding site for optical surveys. 


3. The Simulation of SNe Ia 


The simulation of SN Ia light curves is discussed in 
Section 3.1. We take atmospheric extinction, seeing, telescope 
efficiency, weather, and lunar phase into consideration. We 
explore the impact of exposure time on the accuracy of light 
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Figure 1. The current total efficiency of Mephisto in u, v, g, r, i and z bands. 


measurement and distance measurement in Section 3.2. In 
Section 3.3, we list some well-observed spectra captured 
by LIT. 


3.1. SN Ia Light Curve 


The modified spectral adaptive light curve template (SALT2, 
Guy et al. 2007) model produces a spectral energy distribution 
(SED) for the following simulation. The flux is defined as 
follows 


F(t, A) = xo x [Mo(t, A) + Mit, A)] x exple x CL(A)], 
(2) 


where A and t represent the rest-frame wavelength and time 
respectively, Mo(t, A) and M,(t, A) describe the temporal 
variation of SED, CL()) is the average color correction law, xo, 
x, and c are the fitting parameters of the light curves for SNe Ia, 
xo is the flux normalization parameter, x; is the stretch 
parameter, and c is color. In addition to atmospheric extinction, 
seeing, and telescope efficiency, we consider several other 
factors, including: 

Weather: Lijiang experiences a notable rainy season from 
June to September. The number of nights suitable for 
observation during this period is relatively limited. This paper 
primarily focuses on analyzing the weather patterns during the 
dry season. Based on the statistics of the weather conditions 
(Xin et al. 2020), usually there are about 200 nights in one year 
which are suitable for observations. In this study, we adopt 200 
observational nights per year in the simulation. The occurrence 


of supernovae is uniformly and randomly distributed among 
these 200 days. Furthermore, weather conditions are also taken 
into account. Even during the dry season, there may be 
occasional instances where observation is hindered due to 
factors such as precipitation or high humidity. Therefore, a few 
data points are randomly omitted from the observed data set to 
account for these circumstances. 

Lunar phase: The influence of the lunar phase on observation 
is predominantly concentrated in the u and g bands, which 
correspond to the u, v, and g bands of Mephisto. We 
incorporated the effect of lunar phase into these three bands, 
based on data from Xin et al. (2020), as per the corresponding 
lunar calendar during the time of observation. The observation 
data with and without the influence of the lunar phase are 
presented in Figure 2. As we can see, the influence of the lunar 
phase decreases the number of observable days in the u and v 
bands. 

Data: We generate a batch of SNe Ia with the following 
conditions. In our simulation, the redshifts are less than 0.5. 
The volumetric rate of SNe Ia we adopt in this simulation is 
1.0 x 1074 yr~' Mpe~? according to Rodney et al. (2014). The 
redshift distribution of 4957 SNe Ia with 200 days of 
observation time and a 180 deg? sky area is depicted in the 
upper panel of Figure 5. The absolute magnitude follows a 
normal distribution with a mean of —19.26 and a standard 
deviation of 0.2 (Richardson et al. 2014). The distributions of 
the parameters xı and c obey the following functions (Scolnic 
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Figure 2. This figure simulates an SN Ia with z= 0.1, taking atmospheric 
extinction, seeing, weather, and lunar phase into account. The exposure time is 
130 s. The absolute magnitude in the B-band is —19.3 mag. The top panel 
shows the observation data without the influence of the lunar phase, while the 
bottom panel includes the lunar phase. 
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where A = ./2/7(o, + oy !, and ¥=1.142, 0;=1.652, and 
oo=0.104 for xı, and x = —0.061, 0, = 0.023, o> = 0.083 for c 
(Li et al. 2023). With these distributions, we have compiled a 
comprehensive library of SN Ia light curves, which serves as a 
valuable resource for further analysis. 


3.2. Exposure Time 


We test the photometric accuracy at different exposure times. 
The results are plotted in Figure 3. Due to the sensitivity of the 
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Figure 3. The horizontal axis is exposure time, and the vertical axis is the 
average mag error over the 10 days before and after the B-band maximum. The 
exposure time is from 10 to 200s, with two decimal places retained for the 
result. We select an SN Ia with a redshift of 0.1. 


Table 1 
The Fitting Results Under Different Exposure times, and u Represents the 
Distance Modulus 


Exposure Time H 

20 s 39.02 + 0.21 
60 s 38.84 + 0.11 
80 s 38.83 + 0.10 


130 s 38.82 0.09 


SALT2 fitting to observation points near the maximum 
brightness, we select a total of 10 data points before and after 
the maximum and evaluate the quality of photometric 
measurements based on their errors. Additionally, for bands 
with fewer than 10 data points, we include all available data in 
our analysis. 

As expected, the photometric error generally decreases 
exponentially with increasing exposure time. However, lower 
apparent magnitudes become visible in the u-band with 
increasing exposure time, which leads to an overall increase 
in the average photometric error in the u-band. Only data from 
the g and r bands can be captured when the exposure time is 
limited to 10 s. Data of v, g, r, and i bands are available in all 
four bands when the exposure time is longer than 20s. Data 
from the u band also become observable when the exposure 
time is further increased to 60 s. Data from all six bands can be 
detected when the exposure time exceeds 70 s. An increase of 
10 s in exposure time does not bring significant improvements 
in accuracy at 80 and 130 s. We take into account the SALT2 
fitting accuracy under various exposure times for further 
analysis. Table | displays the results of the fitting. Fitting 
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Table 2 

The Observation Results of SNe Ia with a Redshift of 0.1 under Different Cadences 
SN Ia Ams5(B) Mpg Type(wang) Type(branch) References 
13dy 0.90 —19.65 NV CN Zhai et al. 2016 
09ig 0.90 —19.46 HV nee Marion et al. 2013 
12fr 0.85 —19.49 HV SS Zhang et al. 2014 
15bq 0.82 —19.68 e SS Li et al. 2022 
19ein 1.35 —18.71 HV BL Xi et al. 2022 
l1fe 1.18 —19.40 NV CN Zhang et al. 2016 
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Figure 4. Temporal spectra of different subclasses of SNe Ia at a few days before and after the B-band maximum. The spectrum of 13dy-7 is from Zheng et al. (2013), 


and the others are from LJT. 


results show very small and close errors when the exposure 
time exceeds 60 s. This indicates that the telescope can achieve 
good observation accuracy for SNe Ia with redshifts near 0.1. 
The error in the distance modulus decreases to 0.1 mag when 
the exposure time exceeds 80s. This means the distance 
uncertainty is <4.5% according to the law of propagation of 
uncertainties. We take 80 and 130s into consideration for 
further exploration. 


3.3. Spectral Diversity of Well-observed SNe la 


Spectral follow-up of LJT is crucial for achieving a 
homogeneous SN Ia sample from Mephisto. For example, 
including three spectra at t ~ — 7, 0, + 7 days post-maximum is 
vital in obtaining essential information for diversity studies, 
encompassing the examination of EW, velocity, and velocity 
graduation. Consequently, we can effectively exclude subclass 
samples exhibiting more significant dispersion, thereby 


reducing distance measurement errors arising from scatter 
and enhancing the fitting accuracy of cosmological parameters. 

We present the photometric data of SNe Ia in Table 2, while 
Figure 4 shows spectra of selected well-observed SNe Ia. 
Notably, supernovae with similar Am,; values may exhibit 
varying absolute magnitudes, e.g., SN 2013dy and SN 2009ig. 
Through the measurement of the velocity of Sill A6355, SN 
2013dy is classified into the normal-velocity (NV) type, while 
SN 2009ig is classified into the high-velocity (HV) type in the 
Wang diagram (Wang et al. 2009). The Ams and Mz of SNe 
2012fr and 2015bq are similar. They are both divided into 
shallow silicon (SS) due to the small EWs (<60 A) of Sill 
6355 (Branch et al. 2006). However, SN 2012fr shows 
narrower Sill A6355 and Call IRT profiles, with generally 
deeper absorption lines than the 91T-like events. In the early 
spectra, there are strong high-velocity features of Sill 46355 
that are not seen in other 91T-like events (Zhang et al. 2014). 
SN 2019ein is a 91bg-like supernova. Its peak brightness is 
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Table 3 
The Number of SNe Ia that can be Observed After Screening at Different 
Cadence and Exposure Times 


Cadence Exposure Time Number 
2 days 130 s 1094 
2 days 80s 762 
3 days 130s 979 
3 days 80s 616 
5 days 130s 566 


5 days 80s 312 


Note. The observation time is one year. 


significantly lower than normal SNe Ia. During maximum, the 
91bg-like events show weak or absent Fe II lines and strong 
absorption lines of intermediate elements. These supernovae do 
not follow the WLR. SN 201 1fe is a typical normal SN Ia. We 
also list some SNe Ia with z>0.05, e.g., SN 2022zsy 
(z~ 0.06), SN 2022aadh (z~0.076), and SN 2022adfs 
(z ~ 0.088). 

The collaboration between LJT and Mephisto enables us to 
acquire a substantial quantity of samples featuring spectra at z 
~ 0.1 for subsequent classification based on multidimensional 
spectral information. This information facilitates the elimina- 
tion of subclasses that significantly deviate from the WLR, 
resulting in a more homogeneous sample set. Additionally, 
distinct subclasses of supernovae exhibit varying Phillip 
relationship coefficients, effectively mitigating the dispersion 
observed in the SNe Ia Hubble diagram. 


4. Calibration and Constraints on Cosmological 
Parameters 


To achieve higher accuracy, we conduct a preliminary 
screening on the data before performing the light curve fitting. 
We remove some supernovae with a limited number of 
observed data points. The filtering criteria are as follows: 


(1) signal-to-noise ratio (S/N) greater than 5; 

(2) The sum of all the bands is at least 15 data points; 

(3) Data are available in at least three different bands; 

(4) After applying the aforementioned selection criteria, we 
retained only the supernovae with successful convergence 
in the SALT2 fitting. 


We refer to SNe Ia that have passed the screening as good- 
samples. Table 3 presents the number of supernovae obtained 
under the observation conditions of a 200 day time span and a 
180 deg? sky area. Mephisto has many scientific goals, such as 
the search for extremely metal-poor stars and transient sources. 
In our simulation, we focus on a 180 deg? sky area for our SNe 
Ia survey, which accounts for about 20% of the observing time. 
Our analysis primarily focuses on the cases with cadence of 2 
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days, 3 days, or 5 days, considering exposure times of 80s 
and 130s. 

Various factors can influence the details of an SN Ia 
explosion, including the precise location of ignition, whether it 
is at the center or off-center, and the dynamical burning 
mechanism (e.g., subsonic deflagration or supersonic detona- 
tion). Additionally, the circumstellar environment in which the 
explosion occurs plays a significant role. These factors 
collectively contribute to the observational diversities and 
result in distinct light curve parameters. The dispersion caused 
by these influential factors can substantially impact accurately 
fitting cosmological parameters. Therefore, prior data calibra- 
tion becomes imperative before undertaking cosmological 
parameter fitting procedures. Our study employed the mini- 
mizing x” methodology (Marriner et al. 2011) to produce the 
best fitting result. 


x (a, B) = X imni — wz) 


i=1 
+ax; — bci — Mpin(@P/(o? + oh), (4) 


where Mi, X1; and c; are the best-fit parameters of the SALT2 
model for the i-th SN Ia. Myin is a constant in each redshift bin. 
The number of bins exceeds 30, with each redshift bin having a 
width of 0.01. The definition of Myin is as follows 


Mpin (Zo) = mx — pz) + axı — Ge, (5) 


where m, = —2.5 log xo, a(z) is the distance modulus, and the 
redshift z, represents the central point of each bin. One of the 
key benefits of this approach is its robustness against variations 
in wavelength bands. The difference between m, and mọ is 
about 10 mag. While a cosmological model is employed to 
calculate u(z) in this process, it has been shown in Marriner 
et al. (2011) that the values of a and @ are not affected by the 
choice of cosmological parameters when a large number of bins 
are used. The definition of Gint is as follows 


Tt = Vn, + a? Va + BV. + 2aVin,,x; = 28Vin..c = 2aBV«.c5 


(6) 


where V represents the intrinsic covariance matrix of 
parameters, Vm,, V,,, and V, are the diagonal elements of the 
covariance matrix, and Vinx, Vn„co and Vac are the corresp- 
onding off-diagonal elements. ø? is defined as follows 


oF = Vma + a? Va; + BVa; ae 2A Vmax 
—2BV maici a 2a pVanci (7) 


Similar to the definition of Ti the covariance matrix consists 
of m, x;, and xo for the i-th SN Ia. Restrictions on 
cosmological parameters can be made based on the observed 
distance moduli yz and the theoretically calculated values of pun. 
We consider the flat ACDM model here. We used the MCMC 
(Foreman-Mackey et al. 2013) method for parameter fitting, 
with particular focus on the values of Hp and Qy. The 
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Figure 5. The top panel shows the original number-redshift distribution of 
simulated SNe Ia (initial data, blue histogram) and selected SNe Ia (selected 
data, red step line). The bottom panel displays the moduli-redshift distribution 
of selected SNe Ia (red points). The solid line in the bottom panel represents the 
theoretical prediction. The figure illustrates a scenario where the cadence is 3 
days and the exposure time is 130 s for red points. The observation time is one 
year for both initial data and selected data. 


likelihood function is defined as follows 


£ i i> H; ? Q js 
Hea ys A Hin Ho n)? (a5 
2:20 O p; 


where u; = m, — Mo(zp) + axı, — bci, and fy, is defined as 
follows 


Lin Zi» Ho, Qu) = Slog(d,/Mpc) + 25. (9) 


Considering a flat universe and discarding the radiation term, 
the luminosity distance dz is defined by 


c( +z) f dz! 
Ho 0 Vua +z) + (1 — Qy) 


The absolute peak magnitude of SNe Ia in our simulation 
follows a normal distribution with a mean of —19.26 and a 
deviation of 0.2 (Richardson et al. 2014). This can lead to the 
production of supernovae that are very bright or very dim. We 
removed data that significantly deviated from the Hubble flow. 


d = (10) 


5. Result 


Figure 5 displays the redshift distribution of the screened 
data. The observation strategy utilized a cadence of 2 days and 
an exposure time of 130 s. The top panel displays the 
distribution of the number of SNe Ia with redshifts. The blue 
histogram represents the original redshift distribution of 
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Table 4 
The Observation Results of an SN Ia with a Redshift of 0.1 under Different 
Cadences 
Observation Strategy Ho Êm 
cadence = 2 days, t = 130 s 68.761939 0.44+0.07 
cadence = 2 days, t = 80 s 68.83191 0.454915 
cadence = 3 days, t = 130 s 69.73.2545 0.361508 
cadence = 3 days, t = 80 s 69.7071 53 0.377913 
cadence = 5 days, t = 130 s 70.177 1-03 0.315011 
cadence = 5 days, t= 80 s 70.58/13 0.265915 


simulated SNe Ia, while the red step line depicts the number 
distribution of SNe Ia after the selection process. The bottom 
panel displays the distribution of distance moduli with redshift. 
The solid line represents a flat ACDM model with Qm, = 0.3, 
Qa, =0.7, and Ho = 70 kms! Mpc’, which are used to 
generate the SNe Ia data. The red points in the chart represent 
the distribution of SNe Ia after the selection process. Under this 
observational strategy, the telescope’s observational limit 
extends slightly beyond a redshift of 0.3. For supernovae with 
a redshift around 0.1, there is not only good observational 
accuracy but also no significant selection bias. However, 
selection bias is apparent when the redshift exceeds 0.26. 

The fitting results are shown in Table 4. The obtained fit 
results do not align with the parameters used to generate the 
SNe data, which will be discussed in the next section. Here we 
focus on the errors rather than the fit values. From the results, it 
can be seen that the fit values have the smallest errors for a 
cadence of 2 days and an exposure time of 130 s. However, 
compared to the strategy with a cadence of 3 days and the same 
exposure time, it increases the observation time by 33%, 
resulting in only a 12.5% decrease in the error of 2, and 
almost no change in the error of Hp. We include 88 high-z 
samples at redshifts higher than 0.5 (Riess et al. 2007; Brout 
et al. 2019). The improvement in precision at this level will be 
washed out when incorporating existing data that include high- 
z samples (see Figure 6). Therefore, taking everything into 
consideration, we recommend adopting the strategy of a 
cadence of 3 days and an exposure time of 130 s for the 
observations. 


6. Discussion and Conclusion 


The results obtained by different strategies indicate a 
deviation from the cosmological parameters in generating the 
supernovae. The discrepancy is attributed to multiple factors. 
First, a sample collected at low-z provides limited constraints 
on cosmological parameters. Figure 7 illustrates the Hubble 
diagram at lower redshifts for different values of Qw and Q4. It 
can be seen that at low redshifts, the curves predicted with 
different Qy and Qa almost overlap, indicating that the 
constraints on Qw and Q; are weak by the low-z samples. Even 
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Figure 6. The MCMC fitting results include high-redshift SNe Ia. The left panel is cadence = 2 days, and the right panel is cadence = 3 days. The exposure times are 


both 130 s with a year of observation time. 
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Figure 7. Hubble diagram for different values of cosmological parameters. Red 
points represent simulated SN Ia data, similar to the data in Figure 5. 


a slight deviation can result in significant discrepancies in the 
fitting values, highlighting the sensitivity of the process. 
Another factor to consider is the diversity exhibited by SNe 
Ia, which stems from variations in ignition locations, explosion 
mechanisms, and environments. This study does not incorpo- 
rate spectral simulations, thus making it difficult to ascertain 
the abnormal supernovae within the simulated data. This lack 
of comprehensive classification of SNe Ia may introduce biases 
into the fitting of cosmological parameters. Moreover, it is 
essential to acknowledge that Mephisto filters are non-standard 


and may introduce certain biases during fitting procedures. 
Lastly, while the light curves of SNe Ia align with actual 
observations, any deviations between fitted and actual values 
can be attributed to potential limitations or inaccuracies within 
the utilized cosmological model for generating SNe Ia at 
varying distances; however, these deviations do not signifi- 
cantly impact overall error estimation. 

This study serves as a preliminary investigation of the 
Mephisto survey. Once Mephisto is fully operational, we can 
combine it with LJT to obtain about three spectra during the 
peak phase of well-observed SNe Ia. Currently, LJT can 
identify 30-50 SNe Ia every year. With the addition of 
Mephisto, the sample size of SNe Ia will be significantly 
increased. The fundamental astrophysics of SNe Ia, including 
their progenitor scenarios and explosion mechanisms, remains 
elusive (Liu et al. 2023). A larger sample size of SNe Ia can 
significantly improve our ability to constrain the complex 
interplay between the progenitor channels, explosion models, 
and observations. 

Although Mephisto may not exhibit a numerical advantage 
in the quantity of SN Ia discoveries compared to other large- 
scale survey telescopes, its observational capabilities are 
formidable when detecting SNe Ia at lower redshifts. 
Furthermore, the Mephisto project prides itself on its excep- 
tional color calibration precision and unique ability to conduct 
simultaneous observations across three spectral bands. As a 
result, these features play a significant role in partially 
mitigating systematic errors arising from variations in site 
weather conditions. 


Research in Astronomy and Astrophysics, 24:045017 (9pp), 2024 April 


We conducted a simulation for Mephisto to investigate the 
accuracy of cosmological parameter estimation under various 
observation strategies, with a fixed observation duration of 
180 deg? and 200 days. The optimal observation strategy we 
have determined is an exposure duration of 130 s and a 
frequency of three days. In our simulation, more than 900 SNe 
can be detected in one year of observation. Implementing this 
strategy ensures commendable photometric precision and 
enables us to detect numerous SNe Ia. Our simulation has 
considered the influence of the lunar phase in the light curve 
and assessed the photometric precision of SNe Ia at z ~ 0.1 for 
various exposure durations. The g, of sources near redshift 0.1 
are about 0.1 mag, with distance uncertainty 4.5%. Compared 
to the existing sample, the distance measurement accuracy has 
improved from 7% to 4.5%. By incorporating spectroscopic 
analysis to exclude abnormal SNe Ia, it is possible to achieve a 
better fit of the Phillips relation coefficients and further enhance 
distance measurement precision. We found that beyond an 
exposure time of 60 s, there is a diminishing improvement in 
photometric accuracy with increasing duration. Once Mephisto 
runs officially for survey purposes, the collaborative observa- 
tions with LJT will substantially increase the sample size of 
SNe Ia and facilitate more comprehensive spectroscopic 
classification, thereby mitigating the scatter induced by the 
diversities. 

Mephisto and LJT can build an extensive SN Ia database, a 
valuable resource for further analysis. This database will 
encompass light curves exhibiting diverse shapes and bright- 
ness levels, along with spectra captured near the peak of the 
light curve. This comprehensive data set will offer a valuable 
tool for investigating the heterogeneity of SN Ia events and 
comprehending their utility as cosmological probes. Through 
meticulous examination of these data, we can derive precise 
distance measurements, delve into the nature of dark energy, 
and unveil fresh insights into the origin and evolution of our 
universe. 
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